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SUMMARY  THE  two-parameter  solid  equations  of  state,  the 

MURNAGHAN  EQUATION,  THE  T AIT  EQUATION,  AND  A  GENERALIZED  FORM 
OF  THE  BIRCH  EQUATION  WERE  SELECTED  TO  FIT  THE  STATIC  PRESSURE-] 
VOLUME  ISOTHERMS  OF  SOLID  AND  LIQUID  HIGH  POLYMERS.  THESE 
EQUATIONS  WERE  FOUND  HAVE  A  WIDE  RANGE  OF  APPLICABILITY. 

THEY  CAN  BE  EXTRAPOLATED  TO  VERY  HIGH  FRESSURE  WITH  GOOD 
ACCURACY.  THE  FIRST  AND  SECOND  DERIVATIVES  APPEAR  TO  BE  WELL 
BEHAVED.  THESE  EQUATIONS  CAN  BE  USED  TO  OBTAIN  ANALYTICAL 
EXPRESSIONS  FOR  THE  GRUNE!  SEN  PARAMETER,  CONTAINING  THE 
VOLUME  OR  PRESSURE  EXPLICITLY  AND  TEMPERATURE  IMPLICITLY. 

THEY  YIELD  VALUES  THAT  SEEM  TO  BE  IN  AGREEMENT  WITH  THE  LATEST] 
THEORETICAL  DEDUCTIONS.  WHEN  USED  WITH  THE  KENNEDY  FUSION 
LAW,  THESE  ISOTHERMAL  EQUATIONS  CAN  BE  USED  TO  PREDICT  MELTING] 
TEMPERATURES  UNDER  VERY  HIGH  PRESSURES.  IT  HAS  A^SO  BEEN 
SHOWN  THAT  THESE  EQUATIONS  CAN  BE  USED  TO  PREDICT  THE  SHOCK 
COMPRESSION  DATA,  AND  A  DISCUSSION  OF  THE  SHOCK  COMPRESSION 
DATA  AS  INFLUENCED  BY  THE  HEGREE  OF  CRYSTALLINITY  OF  THE 
POLYMERS  IS  GIVEN.  ALONG  WITH  THE  PVT  PROPERTiES  OF  POLYMERS, 
THE  SPECIFIC  HEATS  (Cp  AND  C  )  OF  POLYETHYLENE,  POLYCHLOROTR I  — 
FLUOROETHYL ENE,  AND  POLYTETR AFLUOROETHYLENE  UNDER  HIGH 
PRESSURES  ARE  EVALUATED. 


KIY  WOftOft  MURNAGHAN  EQUATION,  TAIT  EQUATION,  MODIFIED  BIRCH 
EQUATION,  M I  E- GRUNEI  SEN  EQUATION,  SIMON  EQUATION,  KENNEDY 
FUSION  EQUATION,  HUGONIOT. 
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SUMMARY 


''The  two-parameter  solid  equations  of  state,  the  Murnaghan  Equation,  the  Tait  Equation, 
and  a  generalized  foftn  of  the  Birch  Equation  were  selected  to  fit  the  static  pressure- 
volume  isotherms  of  solid  and  liquid  high  polymers.  These  equations  were  found  to 
have  a  wide  range  of  applicability.  They  can  be  extrapolated  to  very  high  pressure 
with  good  accuracy.  The  first  and  second  derivatives  appear  to  be  well  behaved.  These 
equations  can  be  used  to  obtain  analytical  expressions  for  the  Griineisen  parameter, 
containing  the  volume  or  pressure  explicitly  and  temperature  implicitly.  They  yield 
values  that  eern  to  be  in  agreement  with  the  latest  theoretical  deductions.  When  used 
with  the  Kennedy  fusion  law,  these  isothermal  equations  can  be  used  to  predict  melting 
temperatures  under  very  high  pressures.  It  has  also  been  shown  that  these  equations 
can  be  used  to  predict  the  shock  compression  data,  and  a  discussion  of  the  shock 
compression  data  as  influenced  by  the  degree  of  crystallinity  of  the  polymers  is  given. 
Along  with  the  PVT  properties  of  polymers,  the  specific  heats  (Gp  md  C^)  of  poly¬ 
ethylene,  polychlorotrifluoroethylene,  and  polytetrafluoroethylene  under  high  pressures 


are  evaluated. 


1.  INTRODUCTION 


There  is  considerable  theoretical  and  practical  interest  in  studying  equations  of  state  of 
solids  under  high  pressure.  For  high  polymers,  there  are  only  scattered  data  available  in 
the  literature  including  the  pioneering  work  by  Eridgman.  A  common  method  of 
correlating  the  data  has  been  the  use  of  a  polynomial— -the  specific  volume  as  a 
function  of  ascending  powers  of  pressure.  The  obvious  disadvantaga  of  this  method  is  that 
the  empirical  equation  of  state  is  limited  by  the  range  of  the  experimental  data.  Extrapolation 
beyond  the  range  becomes  highly  uncertain;  however,  extrapolation  is  sometimes  necessary 
for  estimating  the  pressure-volume  relationship  in  the  region  where  no  experimental  data 
exists.  There  is  no  theoretical  approach  which  yields  results  for  high  polymers  which 
compare  in  accuracy  to  the  experimental  data. 

It  has  been  shown  ^  that  the  bulk  properties  of  organic  high  polymers  under  high  pressures 

behave  similarly  to  the  properties  of  other  non-metalb c  and  metallic  solids.  These  properties 

under  high  pressure  can  be  well  represented  by  one  of  the  several  semi-empirical  equations 

of  state  which  can  be  derived  on  the  basis  of  the  theory  of  finite  strains.  Among  these 

(2) 

equations,  the  Mumaghan  and  Bircb.  equations  are  well  known.  Anderson  has  shown  that 

the  Mumaghan  equation  can  be  used  to  estimate  P-v  isotherms  to  very  high  pressures  many 

times  beyond  the  pressure  range  where  the  parameters  were  determined  with  good  accuracy. 

The  simplest  form  of  the  Birch  equation  has  two  fixed  exponents  and  one  adjustable  coefficient. 

A  generalized  and  modified  form  of  the  Birch  equation  based  on  the  theory  of  finite  strains 

containing  two  adjustable  parameters  was  used  in  the  present  paper.  The  Tait  equation, 

also  a  two-parameter  equation  of  state,  has  customarily  been  used  to  correlate  compression 

(3) 

data  for  liquids.  Nanda  and  Simha'  '  have  applied  the  Tait  equation  to  solid  and  liquid  organic 

(4) 

high  polymers.  MacDonald  has  applied  the  Tait  equation  to  other  solids  and  compared  the 
Tait  and  Mumaghan  equations. 
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These  three  equations  have  been  chosen  to  fit  the  volumetric  data  of  several  organic  high 

(5) 

polymers.  The  volumetric  data  used  were  taken  from  the  work  of  Weir'  '  and  Hellwege, 

/c\ 

Knappe,  and  Lehmann'  .  Weir’s  data  were  taken  at  pressures  up  to  10, 000  atm. ;  while 
those  of  Hellwege  et  al.  cover  a  wide  temperature  range  encompassing  both  solid  and  liquid 
regions  of  the  polymers. 

It  was  found  that  these  three  equations  could  be  extrapolated  to  very  high  pressures  with  good 
accuracy.  Values  of  the  first  and  second  derivatives  obtained  from  these  equations  appear  to 
be  smooth,  well  behaved,  and  in  good  agreement  with  the  values  obtained  from  other  methods. 

On  the  basis  of  the  above  results,  four  areas  were  examined: 

a.  The  specific  heati'  Cp  and  Cy  under  high  pressures  for  three  polymers  wore 
evaluated. 

b.  The  melting  phenomena  of  polymers  under  high  pressure  are  examined  on  the  basis 
of  the  Simon  Equation  and  the  Kennedy  fusion  law.  When  combined  with  the  Kennedy 
fusion  law,  these  isothermal  equations  of  state  can  predict  the  melting  temperature 
of  polyethylene  up  to  30, 000  atmospheres  and  agree  to  within  10  percent  with  the 
experimental  data. 

c.  The  Griineisen  parameter  can  be  directly  calculated  from  these  equations  of  state  by 
invoking  one  of  the  well  known  models  of  solids.  Values  thus  obtained  seem  to  be  in 
good  agreement  with  the  latest  theoretical  deductions. 

d.  The  shock  compression  of  polymers  is  shown  to  be  predictable.  The  few  Hugoniots 
found  in  the  literature  are  examined,  and  it  is  found  that  in  all  cases  insufficient 
information  is  given  on  the  physical  properties  of  the  polymers  and  the  initial 
conditions  of  the  experiments,  which  are  important  in  the  calculation  of  Hugoniots. 
The  observed  decrease  in  slop*  of  the  Hugoniots  at  high  shock  pressures  may  be 
explained  by  the  opposite  effects  of  pressure  and  temperature  on  the  degree  of 
crystallinity  of  the  polymers. 

The  Wide  range  of  applicability  of  these  equations  of  state  seems  to  indicate  a  fruitful  area  for 
theoretical  and  experimental  investigation. 
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2.  ISOTHERMAL  EQUATIONS  OF  STATE 


The  Mumaghan,  Tait  and  a  generalized  form  of  the  Birch  equation  have  been  chosen  to  fit  the 
experimental  data.  The  Mumaghan  equation  can  be  derived  from  the  fact  that  the  bulk 
modulus  of  some  solids  is  a  linear  function  of  pressure  along  an  isotherm, 


-v 


B  +  sP 
o 


(1) 


where  v  and  P  indicate  the  specific  volume  and  the  pressure,  respectively;  is  the  initial 
bulk  modulus  at  atmospheric  pressure,  and  s  is  a  dimensionless  coefficient. 


The  Mumaghan  equation  is  obtained  by  integrating  Equation  1  with  the  pressure  given  as  a 

function  of  (v  /v), 
o 


P 


(2) 


(2) 

Anderson  has  shown  that  the  Mumaghan  equation  is  capable  of  predicting  tne  compression 
of  solids  to  very  high  pressures.  He  further  demonstrated  that  the  necessary  parameters 
can  be  obtamed  from  ultrasonic  measurements.  Because  of  the  inherent  inaccuracy  in  the 
direct  measurements  of  specific  volume  and  the  small  change  of  volume  under  pressure,  it 
is  very  difficult  to  obtain  B^  and  B  from  the  slopes  of  the  P-v  isotherm,  especially  when  the 
pressure  interval  is  large. 


In  the  original  Birch  equation,  the  exponents  of  the  v  /v  terms  had  the  fixed  values  of  7/3  and 

o 

5/3.  Because  of  its  greater  flexibility,  a  modified  and  generalized  form  of  the  Birch  Equation 
given  by  Knopoff'  1  is  preferred.  In  this  equation,  the  pressure  is  given  as 
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P  =  3K  /2  (1  +  a  )|(v  /v) 

O  2  1  O 


2  (1  +  a2)/3 


2  (1  +  a  )/3 

(v  /v)  -  1 

0 


(3) 


where  X  =  -v  (dP/3v)_  at  P  =  0,  a„  is  the  coefficient  of  the  quadratic  term  in  tar  f  train 
o  T  2 

energy  expressed  as  a  function  of  strain,  and  P  and  v  have  the  usual  mearings.  If  Equatic.. 

3  is  rearranged,  and  3  K  /2  (1  +  a)  and  3/2  (1  +  a  )  are  considered  as  two  :  ijustable 

o  2  2 

parameters,  the  following  equation  is  obtained 


V/ v 


1/2  +  (P/A'  +  1/4) 


1/2 


-B’ 


(4) 


with  A'  -  3K  /2  (1  +  aj  and  B'  =  3/2  (1  +  aj. 

O  2  2 


The  Tait  Equation  is  purely  empirical.  It  has  been  used  extensively  to  represent  the  P-v 

isotherms  of  liquids.  Recently  the  Tait  Equation  has  been  applied  to  polymers  and  other 
(3  4) 

solids  ’  .  This  equation  can  be  written  as 


v/v  =  1  -  A  In  [(B  +  P)/(B  +  1)] 
o 


(5) 


In  Equation  5,  A  and  B  are  the  two  adjustable  parameters.  The  constant  A  is  often  regarded 

(3) 

as  a  universal  constant.  Nanda  and  Simha  have  taken  this  view  and  apply  the  constant  A 
found  for  hydrocarbons  to  polymers.  In  order  to  obtain  a  good  fit,  both  parameters  have 
been  considered  as  adjustable  in  the  present  work.  It  is  found  that  the  variation  of  A  is 
rather  large  near  the  glassy  transition  temperature  in  the  case  of  polym ethyl  methacrylate. 

An  attempt  has  been  made  to  fit  these  equations  to  the  experimental  data  by  means  of  the 
method  of  least  squares.  The  general  method,  developed  by  Demingv  ,  which  is  applicable 
to  non-linear  equations,  was  not  found  to  be  helpful  in  this  case.  It  seems  to  require  that  the 
initial  estimate  of  the  values  of  the  parameters  be  very  close  to  the  true  values  in  order  that 


o 


the  linearization  of  the  expressions  for  the  errors  be  valid.  Consequently,  an  iteration 

procedure  was  developed.  First,  a  pair  of  initial  values  of  the  parameters  are  obtained 

from  a  pair  of  well  separated  data  points.  Then  the  sum  of  error  squares  is  minimized  by 

trial  and  error  with  respect  to  each  parameter  alternately,  with  the  other  parameter  held 

""6 

constant.  The  minimum  is  of  the  order  of  10  for  Weir's  high  pressure  data.  It  is  in  the 
order  of  10  for  the  data  of  Hellwege  et  al.  The  sum  of  error  squares  becomes  larger  as 
the  melting  point  is  approached.  In  the  case  of  polyethylene  in  the  melting  region  (where  the 
volume  change  is  rather  sudden),  fitting  of  the  data  could  only  be  done  by  obtaining  separate 
parameters  for  the  liquid  and  solid  phases. 

Tables  1,  2,  and  3  give  the  parameters  in  the  Mumaghan,  modified  Birch,  and  the 
Tait  equations.  Temperature  dependence  of  the  parameters  in  each  of  the  equations  in  the 
case  of  polymethyl  methacrylate  are  plotted  in  Figures  1,  2,  and  3.  The  parameter 

in  the  Murnaghan  Equation,  the  parameter  in  the  modified  Birch  Equation,  and  B  in 
the  Tait  Equation  are  always  decreasing  with  increas.ng  temperature  as  expected.  Both 
coefficients  s  (Murnaghan  Equation)  and  a  (Birch  Equation),  which  are  related  to  the 
compressibility,  remain  fairly  constant  at  temperatures  far  below  the  transition  temperature 
and  then  rise  and  f:.ll  precipitously  on  going  through  the  transition. 

The  parameter  s  in  the  Murnaghan  equation  for  solid  and  liquid  atactic  (amorphous)  polystyrene 
is  plotted  in  Figure  4.  The  shape  of  the  curve  for  the  parameter  s  is  similar  to  that  of  the 

(9) 

specific  heat  data  of  Karasz,  Bair,  and  O'Reilly  .  The  cause  of  the  abrupt  rise  of  s  near 
o  (3  \ 

250  C  is  not  clear.  Nanda  and  Simha  '  suspect  degradation  or  loss  of  material  at  this  high 
temperature. 

Any  of  the  three  equations  of  state  can  be  used  to  calculate  the  specific  volumes  of  polymers 
with  higher  accuracy  than  can  Weir's  polynomials,  even  in  the  range  where  the  polynomials 
were  fitted  to  accurately  represent  the  experimental  data.  In  Table  4,  the  specific  volumes 
of  polytetrafluoroethylene  at  30°C  and  at  pressures  from  1000  to  10,000  atmospheres  are 
compared  with  values  calculated  from  the  Murnaghan,  Birch,  Tait,  and  the  polynomial 


Table  1.  Parameters  in  the  Mumaghan  Equation  of  State  for  High  Polymers 


Polymer 

B 

0 

(atm) 

s 

V 

0 

T  (°K) 

Static 
Pressure 
Range  (atm) 

Low  Density 
Polyethylene 

30063.2 

25138.4 

9.4902 

8.4305 

1.0895 

1.1269 

293.2 

343.2 

10,000 

10,000 

High  Density 
Polyethylene 

(a) 

40148.4!  : 
35034. 6 (a) 

16.9992 

12.2738 

1.028 

1.042 

292.6 

332.8 

Poly  methyl 
Methacrylate 

39712. 3  j3; 
28191. 3(a) 

9.0004 

10.9414 

0.842 

0.862 

293.2 

373.2 

Polystyrene 

36388. 4 fj 
14073.2 

8.5179 

15.7109 

0.955 

0.994 

293.2 

399.2 

2,000^1 

2, 000  a 

Potychloro- 
trifl.  oroethylene 

39146.6 

33993. 5 

10.4331 

9.  8502 

0.4662 

0.4745 

293.2 

343.2 

10,000 

10,000 

Polytetrafluoro- 

ethylene 

29267 

21040 

14.9218 

12.8603 

0.4484 

0.4484 

293.2 

293.2 

1  to  5,000 
6,000  to  10,000 

Polyvinyl 

Fluoride 

416S7.7 

31102.5 

7.8864 

9.3506 

0.7956 

0.8135 

293.2 

343.2 

10,000 

10,000 

Polyvinylidene 

Fluoride 

34757.8 

32208.2 

9.5880 

9.3077 

0.6439 

0.6582 

293.2 

343.2 

10,000 

10,000 

Polyviny 

Alcohol 

58233, 0 

46258.8 

10.7629 

11.0647 

0.7720 

0.7861 

293.  2 

353.2 

10,000 

10,000 

Nylon  6-6* 

Nylon  6-10 

56688. 2, a} 
55470. 7(a) 

7.7507 

7.8852 

0.858(b) 

0.917 

m 

(a) 

40,000  ' 
40, 000  a 

(a)  the  unit  is  in  Kg/cm 

(b)  Based  on  average  density  of  amorphous  and  crystalline  samples. 
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Table  2.  Parameter  in  the  Modified  Birch  Equation  of  State  for  High  Polymers 


Polymer 

Ko 

(atm) 

a2 

Temperature 

(°K) 

Low  Density 

27861.4 

4.77267 

293.2 

Polyethylene 

22116.0 

3.83297 

353.2 

High  Density 

42770. 0(a) 

7.0098 

292.6 

Polyethylene 

9063.7 

10.1728 

403.3 

Polymethyl 

39813.  9(a) 

3.54545 

293.2 

Methacrylate 

27704.4 

5.08729 

373.2 

36284.  7 (a) 

3.4300C 

293.2 

Polystyrene 

17976.6 

8.9575 

373.8 

Polychloro-tri- 

39353.3 

4.5555S 

293.2 

fluoroethylene 

31634.9 

4.97655 

353.2 

Polytetrafluoro- 

27741.8 

7.7944 

293.2  (0  -  5000  atm) 

ethylene 

.  . . - 

18801.6 

7.0129 

293.2  (5000  -  10000  atm) 

Polyvinyl 

39714.1 

3.64910 

293.2 

Fluoride 

29396.5 

4.54066 

343.2 

Polyvinyl 

57311.2 

4.  93054 

203.2 

Alcohol 

-  j 

49088.0 

4.25164 

353.2 

Nylon  6-6  (oriented) 

(a) 

52578.3;  : 

3.69220 

300 

Nylon  6-6  (unoriented) 

62392.  7  jaj 

3.56303 

Nylon  6-10 

50865.  7'a 

3.83297 

300 

Table  3.  Parameters  in  the  Tait  Equation  of  State  for  High  Polymers 


Polymer 

A 

B  (atm) 

Temperature  (°K) 

Low  Density 

0.086812 

2510.1 

293.2 

Polyethylene 

0.097213 

2383.0 

343.2 

High  Density 

0.055221 

2310. 9(a) 

292.6 

Polyethylene 

0.  054046 

1809.7 

332.8 

Polymethyl 

0.098701 

3911.? (a) 

293.2 

Methacrylate 

0.  080562 

2253.7 

373.2 

Polystyrene 

0.  090000 

(a' 

3200. 0V  ' 

293.2 

Polychloro- 

0.086480 

3430.0 

293.2 

trifluoroethylene 

0.084714 

2777.4 

353.? 

Polyt  et  rafluoro  - 

0.0558? 

1523.5 

295.2  (1  to  5000  atm) 

ethylene 

0.06367 

1216.8 

293.2  (6,  000  to  10,  000  atm) 

Polyvinyl 

C. 19650 

4400. 0 

293.2 

Fluoride 

0. 088401 

2656.8 

353.2 

Polyvinylidene 

0.088353 

3019.4 

293.2 

Fluoride 

0.090014 

2834.2 

343.2 

Polyvinyl 

0.081121 

4676.8 

293.2 

Alcohol 

0.077595 

3513.8 

353.2 

(a) 


The  unit  of  B  is  in  Kg/cm 
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Polymethyl  Methacrylate 
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Figure  3.  Parameters  A  (Left  Scale)  and  B  (Right  Scale)  in  the  Tait  Equation  of  State  for  Polymethyl 

Methacrylate 
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Figure  4.  Parameters  BQ  (Left  Scale)  and  s  (Right  Scale)  in  the  Murnaghan  Equation  of  State  for  Solid  and 

Liquid  Polystyrene 


Table  4.  Comparison  of  Specific  Volumes  (cm3/g)  of  Polytetrafluoroethylene  at  30°C  with 
the  Calculated  Values  by  Polynomlnal,  Murnaghan,  Birch,  and  Talt  Equations  of  State 


P,  atm 

Expt'l 

1000 

0.4367 

2000 

0.4289 

3000 

0.4228 

4000 

0.4174 

5000 

0.4120 

6000 

0.3985 

7000 

0.3942 

8000 

0.3907 

9000 

0.3879 

10000 

0.3856 

Mumaghan  |  Birch 


0.4397 

0.4308 

0.4219 

0.4181 

0.4134 

0.3981 

0.3921 

0.3905 

0.3876 

0.3855 


0.4376 

0.4286 

0.4221 

0.4170 

0.4128 

0.3983 

0.3944 

0.3910 

0.3880 

0.3853 


0.4373 

0.4284 

0.4221 

0.4171 

0.4130 

0.3982 

0.3944 

0.3910 

0.3880 

0.*?53 


Talt 


0.4374 
0.4285 
0  4220 
0.4170 
0.4129 
0.3981 
0.3944 
0.3910 
0.3880 
0. 3853 


equations  of  state.  In  the  transition  region  (1000  to  5000  atmospheres),  all  the  three 
isothermal  equations  give  more  accurate  predictions  than  does  Weir's  12-parameter 
polynomial.  Even  in  the  range  6000  to  10, 000  atmospheres  where  the  polynomial  fits  the 
data  very  well,  these  equations  are  still  more  accurate. 

2 

The  compression  data  (v/v  )  of  polyethylene  and  polystyrene  up  to  40,000  Kg/cm  ,  published 

dO)  ° 

by  Bridgman'  are  compared  in  Table  5  with  values  calculated  from  the  generalized 
Birch  equation.  The  parameters  used  in  the  calculations  were  evaluated  by  fitting  the  data 

O 

of  Hellwege  et  al.  up  to  only  2000  Kg/cm“.  The  agreement  might  have  been  better,  if  the 

o 

"room"  temperature  of  Bridgman's  experiment  had  been  specified.  In  this  calculation  20  C 
has  been  used. 

The  derivatives  of  these  equations  yield  smooth  and  reasonable  values.  The  compressibility 
of  polyethylene  at  20°C  is  calculated  by  means  of  the  Murnaghan  equation  and  compared  to 
that  obtained  from  graphical  differentiation  of  the  P-v  isotherm  and  with  the  compressibility 
calculated  from  the  polynomial  equation.  This  comparison  is  given  in  Figure  5. 
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Table  5.  Comparison  of  Calculated  (by  Birch  Eq. )  v/v  of  Polyethylene,  Polystyrene 
up  to  40,000  Kg/'cm^  with  Bridgman’s  Experimental'*®)  Results 


Pressure 

Polyethylene 

_ 

Polystyrene 

P 

T77  1 

% 

2 

o 

° 

(Kg/cm  ) 

(calc.) 

(Expt’l) 

Dev. 

Dev. 

2500 

0.9378 

0.  9447 

0. 7? 

0.9473 

0. S591 

1.23 

5000 

0.9025 

0.9126 

1.11 

0.9125 

0.  92".  3 

1.63 

10000 

0. 8582 

0.8681 

1.14 

0. 8656 

0.  8800 

1.64 

15000 

0.8293 

0. 8379 

1.03 

0. 8334 

0.  8481 

1.73 

20000 

0.8078 

0.8148 

0.  86 

0.8090 

0.  3244 

1.87 

25000 

0.7907 

0.7960 

0.66 

0.7893 

0.  8054 

2.00 

30000 

0. 7766 

0. 7306 

0.51 

0.7728 

0.  7894 

2.10 

35000 

0.7646 

0.7679 

0.44 

0.7587 

0.7760 

2.23 

40000 

0.7541 

0.7570 

.  . 

0.39 

0.7464 

0.7641 

2.32 
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Figure  5.  Comparison  of  Calculated  Compressibility  of  Polyethylene  with  Values  Obtained  from  Gr  aphical  Slope 

and  Values  Calculated  from  Polynomial  Empirical  Equation  of  State 


3.  SPECIFIC  HEATS  OF  HIGH  POLYMERS  UNDER  HIGH  PRESSURES 


The  thermodynamic  expression  for  the  variation  of  Cp  with  pressure  and  the  relationship 
between  C?  and  Cy  are  well  knovm.  They  are  given  as, 

fP  2  > 

Cp  (P)  =  0^  (P  =  O’  -  T  I  (5  v/8  T")  dP  ( 

•'o 


Cp  (P)  -  Cv  (P)  =  Tvff/D 


where  a,  and  /l  are  the  thermal  expansion  coefficient  and  the  compressibility,  respectively. 
The  quantitives  P,  T,  Cp,  C^,  and  v  have  the  usual  meanings.  The  thermal  expansion 
coefficient  a  and  the  second  derivative  of  volume  with  respect  to  T  for  branched 
polyethylene,  polytetrafluoroethylene,  and  polychlorotrifluoroethylene  can  be  obtained 
as  a  function  of  P  from  the  polynominal  equations  by  Weir^. 


The  volumetric  data  reported  by  Hellwege  et  al.  have  only  three  significant  figures. 

,,  -4 

Since  the  thermal  expanison  coefficient  is  in  the  order  of  10  ,  it  is  difficult  to  obtain 

2  2 

reasonable  values  of  This  quantity  Oln  v/}  T)p  and  of  (}  v/}  T  )p.  At  the  present 
time,  the  specific  heat  of  polystyrene  and  polymethyl  methacrylate  under  high  pressure 
cannot  be  obtained  from  equations  (6)  and  (7). 


The  specific  heats  of  many  high  polymers  have  been  measured  under  atmospheric 
pressure.  In  Table  6  there  are  assembled  the  data  sources,  types  of  material,  and 
the  range  of  measurements  of  several  polymers  which  are  of  the  most  interest  in  this 
work. 
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Table  6.  Data  Sources  of  Specific  Heat  (Cp)  of  High  Polymers  under  Atmospheric  Pressure 


Polymer 

Type  of 

Polymer 

Temperature 

Range 

n 

Polyethylene 

Linear  and  branched 

0  o 

90  to  433  K 

(11) 

Polyethylene 

Annealed,  drawn, 
granular  and  sheet 

-20°  to  200°C 

(12) 

Polyethylene 

Low  pressure 

53°  to  145°C 

(13) 

Polychloro- 

trifluoroethylene 

Quenched  and 
slow  cooled 

0°  to  240.  8°C 

(14) 

Polytetrafluoro- 

Powdered  and  quenched 

o  o 

346  to  725  K 
o  o 

(15) 

ethylene 

Powered  and  drawn 

-25  to  120  C 

(16) 

Polytetrafluoro- 

ethvlene 

Commercial 
sheet  sample 

-130°  to  80°C 

(17) 

Polvmethyl 

Methacrylate 

(Data  by  Sochava,  et  al) 

10°  to  300°K 

(18) 

Polyvinyl  Alcohol 

(Data  by  Sochava,  et  al) 

10°  to  245°K 

(19) 

Polystyrene 

Isotactic 

o  0 

0  to  310  K 

(20) 

Polystyrene 

Atactic  -ann  ea1  ed , 

293°  to  378. 61°K 

(9) 

as  received 

89.43°  to  473. 14°K 

(9) 

Isotactic -amorphous, 

301. 9°  to  521. 4°K 

(9) 

annealed 

305°  to  526°K 

(9) 

Polystyrene 

Atactic  and  isotactic 

-50°  to  280°C 

(21) 

Polyvinyl  Chloride 

With  10,20,30r:  plasticizer 

-30°  to  80°C 

(17) 

Polyvinyl  Chloride 

Powdered  GE  L-38 

-20°  to  120°C 

(22) 
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The  specific  heats  Cp  and  Cv  under  pressures  up  to  10,  000  atmospheres,  have  been 
obtained  for  polyethylene,  polychlorotrifluoroethylene,  and  polytetrafluoroethylene 
from  Weir's  polynomials.  They  are  presented  graphically  in  Figures  6  through  12. 

In  Figures  6  and  7,  the  specific  heat  of  polyethylene  is  plotted  as  a  function  of 
temperature  and  as  a  function  of  pressure.  Some  errors  in  the  values  of  the  specific 
volumes  of  polyethylene  were  found  in  Weir’s  data  between  1  and  2000  atmospheres 
and  above  60°C.  These  errors  have  been  corrected  by  comparison  with  the  data  of 
Hellwege  ec  al.  The  material  used  by  Weir  was  a  bracked  type  polyethylene  which  has 
a  melting  point  of  about  110°C.  The  atmospheric  specific  heat  Cp  of  this  material  has 
been  measured  by  Dole  et  al^^.  This  data  has  been  used  in  the  Cp  (0)  in  Equation  (7) 
to  obtain  the  Cp's  shown  in  Figures  6  and  7. 

/P  2  ,  2 

The  effect  of  pressure  on  C_,  i.e. ,  T  /  (8  v/8  T  )p  dP  for  polytetrafluoroethylene 

P  Jo  -4  _3 

and  polychlorotrifluoroethylene  are  in  the  order  of  10  and  10  ,  respectively. 

Compared  to  the  magnitude  of  the  specific  heat  itself,  this  effect  is  negligible. 

The  quantity  (Cp  -  C^)  of  polytetrafluoroethylene  is  plotted  as  a  function  of  temperature 
in  Figure  10.  Two  first  order  transitions  have  been  reported  between  10  and  40°C 
and  1  and  6,000  atmospheres.  This  transition  is  reflected  in  the  (Cp  -  C^)  values. 

Above  6,000  atmospheres  (Cp  -  C^)  becomes  very  small. 

The  quantity  (Cp  -  C^)  for  polychlorotrifluoroethylene  (Kel-F)  as  a  function  of  temperature 
and  of  pressure  is  plotted  in  Figures  11  and  12,  respectively.  Figure  11  shows  that 
(Cp  -  C^)  is  linear  in  temperature  up  to  9,000  atmospheres.  In  Figure  12,  the  three 
isotherms  cross  each  other  twice,  once  between  3,000  to  4,000  atmospheres  and  once 

3 

near  8  x  10  atmospheres.  The  one  atmosphere  curve  does  not  show  an  effect  of  the 
glass  transition  at  52°C  and  one  atmosphere  which  has  been  confirmed  by  Hoffman  and 
Weeks' (23). 
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Figure  6.  Pressure  Dependence  of  the  Specific  Heat,  Cp  of  Polyethylene  (low  density)  at  Constant 

Temperatures 
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Figure  s.  Temperature  Dependence  of  (C  -  C  )  of  Polyethylene  Along  Isobars 
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Figure  9.  Pressure  Dependence  of  (C  -  Cv)  of  Polyethylene  Along  Isotherms 
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Figure  10.  Temperature  Dependence  of  (Cp  -  C  )  of  Polytetrafluoroethylenc  Along  Isobars 
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4.  THE  MIE-GRUNEISEN  EQUATION  OF  STATE 


The  thermal  energy  of  crystalline  materials  consists  of  the  vibrational  energy  of  atoms  in 
the  lattice.  The  Gruneisen  parameter  car.  be  expressed  in  terms  of  specific  volume  v  and 
the  freouencv  • .  of  the  ha r monk,  oscillators, 

i 

7  -  -  (d  In  v  /d  In  v)  (3) 

av. 

The  Mie-Griineisen  equation  oi  state  may  be  written  as: 

P  -  P  =  7  (E  -  E  )  v  (9) 

o  o 

where  P  and  E  are  pressure  and  energy  of  a  reference  state,  and  r  is  the  Griineisen 
o  o 

parameter. 

Differentiating  P  in  Equation  9  with  respect  to  E  at  constant  volume  and  substituting 
dE  =  dT.  7  can  be  expressed  in  terms  of  experimental  quantities  as: 

7  =  vof/C  p  GO) 

v 

where  a  and  p  are  the  thermal  expansion  coefficient  and  the  compressibility  respective] y. 

From  equations  of  state  and  the  exDerirnental  data  discussed  above,  7  car.  be  easily  calculated. 
The  values  of  7  for  several  polymers  have  been  evaluated  in  Reference  1.  Figures  13  and 
14  present  the  Gruneisen  parameter  of  Kel-F  as  functions  of  temperature  and  volume, 
respectively.  The  volume  dependence  of  this  parameter  is  well  known,  but  the  temperature 
dependence  has  not  .  een  clearly  established.  In  the  case  of  Kel-F  in  this  limited  temperature 
range,  the  Gruneisen  parameter  decreases  significantly  with  temperature  increase.  Recently, 
Barker'  has  concluded  that  the  Gruneisen  parameter  of  polymers  is  ««  4  and  the  77 /dT  >  0. 
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Based  on  a  simple  model.  Slater 


(25) 


has  given  an  expression  for  the  Gruneisen  parameter  as, 


r  =  -2/3  -  v  02  PA  v2)t/2  (b  P/S  V)T 


(11) 


This  expression  employs  only  the  isothermal  P-v  property  of  a  solid.  This  formula,  as 

2/3 

modified  by  Dougdale  and  McDonald,  contains  derivatives  of  Pv  .  It  has  recently  been 

/OfJ  \ 

generalized  by  Al'tshuler  ,  who  gives 


=  -(2/3-z/S)  -  v  (S2  Pv2z/3/Sv2)/2(SPv2z/3/Sv)r 


(12) 


which  includes  the  above  two  as  special  cases.  The  Slater  theory  corresponds  to  z  =  0,  and 
the  Dougdale-McDonald  theory  corresponds  to  z  --  1. 

These  formulas,  which  do  not  require  calorimetric  measurements,  offer  convenient  means 
to  obtain  the  Gruneisen  parameter  from  any  one  of  the  three  isothermal  equations  of  state  that 
have  been  fitted  to  the  experimental  data  in  this  work.  Based  on  Slater’s  theory,  the 
expressions  for  the  Gruneisen  parameter  from  each  of  the  three  equations  of  state  are  given 
as  follows: 

Mumaghan  Equation 


=  (s  -  l)/2  -  2/3 


(13) 


Modified  Birch  Equation 


-2/3  -  (1/2B')  [{B’  +  2)  (v/v 


,-l/B' 


(B'  +  l)/2 


C(v/v  ) 


-1/B' 


-  1/2] 


(14) 


or 
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r  =  -2/3  +  (1/4B*)  (P/A1  +  1/4) 


+  (B'  +  2)/2B* 


(15> 


-1/2 


Tait  Equation 

F  =  (v/v  )/2A  -  2/3  (16) 

o 

The  meanings  of  the  parameters  s,  A',  B'  and  A  have  been  defined  in  Equation  2,  Equation 
4  and  Equation  6. 

In  Equation  13,  T  is  dependent  only  on  s,  which  increases  with  increasing  temperature  in 
the  cases  of  polymethyl  methacrylate  and  polystyrene  (Figure  1  and  Figure  4).  At  lower 
temperatures,  remote  from  the  melting  zone,  F  is  in  the  order  of  4. 

In  Equations  14  and  15,  the  dependence  of  F  on  volume  or  pressure  is  explicit,  while  the 
dependence  on  temperature  is  through  A'  and  B'.  Since  the  first  factor  (2B1)  1  is  increasing 
with  temperature  (Figure  2),  examination  of  the  relative  magnitudes  of  the  terms  in  this 
equation  shows  that  the  derivative  dF/dT  is  positive  for  polymethylmethacrylate. 

From  Equation  16,  the  temperature  derivative  of  the  Griineisen  parameter  is: 

dT/dT  =  -  (v/v  )  (dA/dT)/2A2  -  [dv/dT  -  (v/v  )  (dv  /dT)  ]  (2Av  f 1  (17) 

O  0  0° 

If  the  Tait  parameter  A  is  independent  of  temperature,  dT/dT  >  0  at  high  pressures  because 
dv  (P)/dT  <  dv^  (P  =  0)/dT.  If  v  is  decreasing  with  increasing  temperature,  as  in  the  case 
of  polymethylmethacrylate  above  50°C  (Figure  3),  the  first  factor  on  the  right-hand  side  of 
Equation  17  is  also  positive.  If  A  %  0. 1  and  dA/dT [dv/dT  -  (v/v0)  (dv  /dT)],  the  first 
term  on  the  right-hand  side  of  Equation  4-10  will  be  greater  than  the  second  term. 
Consequently,  dT/dT  >  0. 
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5.  MELTING  OF  ORGANIC  HIGH  POLYMERS  AT  HIGH  PRESSURES 


The  Simon  equation  has  been  used  extensively  to  represent  the  melting  temperature  of 
solids  as  a  function  of  melting  pressure: 


P 

m 


,  o  c 
(T  /T  ) 

'  m  nr 


(18) 


where  P^  (atmosphere)  and  T^  (°K)  are  the  melting  pressure  and  temperature  respectively. 
If  P^  =  1  is  taken  as  reference  pressure,  then  T^  is  the  melting  point  of  the  solid  at  1 
atmosphere. 


The  melting  temperature  of  polymers  increase  with  pressure  as  is  the  case  of  other  solids  . 
On  the  other  hand,  the  melting  zone  of  high  polymers  covers  a  rather  large  temperature 
range.  The  amorphous  constituent  softens  gradually,  and  the  crystallites  of  increasing  size 
melt  at  progressively  higher  temperatures. 

The  melting  transition  zone  of  polystyrene  (atactic,  ai  -phous) ,  polyvinylchloride, 

polymethyl  methacrylate,  and  two  types  of  polyethylene  (high  and  low  densities)  has  been 

investigated  by  Hellwege  et  al.  ^  under  pressures  up  to  2000  Kg/cm2.  The  melting  of 

several  types  of  polyethylene  under  pressures  up  to  4,  000  atmospheres  has  been  investigated 
(27) 

oy  Matsuoka  .  The  melting  of  a  type  of  high  density  polyethylene  har  been  studied  by 

Osugi  and  Hara^  un  to  30.  000  atmospheres  where  the  melting  temperature  is  400°C. 

(29) 

McGeer  and  Duns'  ;  obtained  data  on  the  melting  of  Teflon  under  pressures  up  tc  815 
atmospheres. 


Based  on  these  results,  the  parameters  for  polymers  in  the  Simon's  equation  have  been 
evaluated.  These  parameters  are  presented  in  Table  7.  A  list  of  the  Simon  constants  for 
solid  elements  and  compounds  has  been  given  by  Babb^°\  The  constants  for  polymers 
presented  in  Table  7  are  comparable  with  those  'or  some  of  the  solids  studied  by  Babb. 


Table  7.  Parameteis  of  the  Simon  Equation  for  Polymers 


Polymer 

a 

(atm) 

c 

Pressure  Range 
(atm) 

Teflon 

1086 

3.329 

1  to  615 

Polystyrene 

19961 

0.5928 

1  to  2000 

Low  Density 
Polyethylene 

4369 

3.261 

1  to  4000 

High  Density 
Polyethylene 
(Marlex  9) 

6742 

2.115 

1  to  2000 

High  Density 
Polyethylene 
(Marlex  50) 

3100 

4.5 

1  to  30,000* 

Polyvinyl  Chloride 
(Hart) 

18869 

1.27 

1  to  2000 

Polymethyl 

Methacrylate 

7526 

2.083 

1  to  1500 

*  J.  Osugi  and  K.  Hara  Rev.  Phvs.  Chem. (Japan)  36,  28,  1966 


Salter  has  pointed  out  that  the  Simon  constant  c  is  related  to  the  Gruneisen  parameter  by 
the  following  equation: 


,(41) 

c  =  (r  +  1/6)  (T  -  1/3)  (13) 

which  would  enable  theoretical  calculation  of  c  if  the  Grh'neisen  parameter  is  known  more 
accurately. 

'31) 

Recently,  Kennedy'  has  discovered  empirically  that  the  melting  temperature  is  a 
linear  function  of  the  isothermal  compression  (Av/v^)  at  normal  temperature: 


t  (P)  =  t°  (1  +  kAv/v  ),  °C  (20) 

m  m  o 

In  subsequent  articles^32  Kennedy  et  al.  have  shown  that  for  the  majority  of  solids 

Equation  20  holds  very  well.  For  a  few  solids,  such  as  solid  argon  and  methane,  the 
deviations  observed  might  be  attributable  to  quantum  effects  or  the  uncertainties  in  the 
compression. 

Gilvarry^^  and  Vaidya  and  Gopal^^  have  derived  the  relations  between  the  constant  k  in 
Equation  20  and  the  Gruneisen  parameter  T,  and  between  k  and  the  constant  c  in  Simon 
Equation: 

k  =  2  (T  -  1/3)  -  (c  -  l)"1  (21) 

The  available  melting  data  of  polymers  have  been  assembled  to  test  the  Kennedy's  fusion 
equation.  Within  the  range  of  the  available  data,  Equation  20  seems  to  be  applicable  to 
polymers.  The  results  for  six  polymers  are  plotted  in  Figure  15,  and  the  constant  k  is 
presented  in  T;ible  8.  The  compression  data  for  polytetrafk’oroethylene  and  low  density 
polyethylene  were  based  on  the  data  by  Weir^.  Both  the  melting  and  the  compression 
data  of  polystyrene,  polyvinyl  chloride,  and  poiv^thyl  methacrylate  were  based  on  the 
data  of  Hellwege  et  al.  ^ . 
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Figure  15.  Application  of  the  Kennedy  Fusion  Law  to  Organic  High  Polymers 


1  thle  8.  Constanta  in  the  Kennedy  Fusion  Equation 


1 

Polymer 

k 

o 

t  <  C> 
m 

Maximum 

4v/v 

n 

Polytetr  afluo  r  o - 
ethylene 

26.78 

325 

0.0108 

Low  Density 
Polyethylene 

9.27 

111. 

0.0829 

High  Density 
Polyethylene 
(Marlex  9) 

10.33 

136 

0.  0276 

Polystyrene 

19.  75 

1 9. 6* 

0.  0441 

Polyvinyl 

Chloride  (Hart) 

11.02 

71.5 

0.  0395 

Polymethyl 

Methacrylate 

11.02 

i 

i 

i 

— — - -  _ — , 

108.1 

! 

0.  0319 

*Value  obtained  by  extrapolation  to  1  atmosphere 


By  combining  Equation  20  with  the  isothermal  equations  of  state,  the  melting  temperature 

can  be  directly  predicted.  The  melting  temperature  of  high  density  polyethylene  (Marlex  50) 

measured  by  Osugi  and  Hara^  '  at  pressures  up  to  30, 000  atmospheres  is  compared  with 

the  calculated  melting  temperature  by  means  of  Equation  20  in  Table  9.  The  compression 

Av/v^  was  calculated  from  the  isothermal  equations  of  state  Equations  2,  3,  and  6. 

The  constant  k  in  Table  8,  and  the  constants  in  these  equations  were  determined  for 

2 

pressures  up  to  2000  Kg/cm  for  Marlex  9  which  is  different  from  Marlex  50  in  the  degree 
of  crystallinity  and  the  number  of  methyl  groups  (CH  )  per  100  carbon  atoms.  Below  10,  000 

O 

atmospheres,  all  those  equations  yield  values  agreeing  with  the  experimental  data  within  1 
percent.  Above  10,  COO  atmospheres  the  difference  between  the  calculated  and  experimental 
melting  temperature  is  greater.  Both  the  Birch  and  the  Tait  equations  give  values  which 
deviate  from  the  experimental  by  less  than  10  percent  at  30,  000  atmospheres.  The  value 
from  the  Murnaghan  equation  differs  slightly  more.  10.2  percent,  at  30,000  atmospheres. 
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Table  9.  Comparison  of  Measured  and  Predicted  Melting  Temperature 
of  High  Density  P  ylene  by  Means  of  Isothermal  Equations  of 
St?1-  md  Kennedy  Fusion  Law 


Pm,  Atm 

V  °K 

(Measured) 

Predicted  T  ,  °K 
m 

Birch  Eq 

Murnaghan  Eq. 

T^t  Eq 

1 

409 

409 

409 

409 

1,000 

431 

437 

437 

437 

3,000 

468 

474 

472 

473 

5, 000 

509 

501 

494 

498 

10, 000 

551 

545 

529 

539  ! 

20,000 

630 

597 

567 

1 

585 

30,000 

671 

- — - - IL 

630 

590 

613 
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0.  SHOCK  COMPRESSION  OF  HIGH  POLYMERS 


Shock  compression  is  a  means  of  obtaining  semi-  empirical  equations  of  state  for  matter 
under  high  pressures  and  high  temperatures.  Properties  of  maxter  under  such  conditions 
are  determined  by  the  interactions  and  thermal  vibrations  of  its  component  molecules,  atoms, 
ions,  and  the  excitation  of  the  electrons.  Static  compression  under  isothermal  conditions 
serves  to  investigate  the  forces  of  interaction,  and  the  tbtfmal  energy  nt  the  prevailing 
temperature.  For  a  given  initial  state,  the  isothermal  pressure  can  be  calculated  for  a 
specified  volume,  v,  by  use  of  the  isothermal  equations  of  state  developed  above.  If  the 
pressure  due  to  temperature  rise  can  be  determined,  then  the  prediction  of  the  shock 
comr  rssion  curve,  better  known  as  Hugouiot,  is  possible. 

In  a  general  r.on-iscthermal  compression, 

dP  -  £P  i  v)T  dv  +  dT  (22) 

with  use  of  the  well  known  relations 

«=P.  =T'  -  <5  In  v  cT)p  (b  In  v 

Equation  6-i  can  be  rearranged  to  give 

-dP.'v  (cP.bv)  =  -  dv  v  r  In  v-sT)  dT  (23) 

Integration  of  Equation  23  will  be  easy  if  the  thermal  expansion  coefficient  o;  =  (c  In  v/cT) 
is  assumed  constant,  ix  has  been  shown'1  \  based  on  such  an  assumption  of  constant  o  ,  that 
a  non-isothermal  Murnaghan  Equation  can  Ire  derived, 
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(24) 


P  =  (B  /s)  [(v  /v)  -  l]  +  B  (v  /v)  a  (T  -  T  ) 

0  0  O  0  0 


where  and  s  are  the  parameters  in  the  Mumaghan  Equation  of  state  and  the  other  symbols 
have  the  usual  meanings.  By  equating  the  Rankin e- Hu goniot  Equation  to  the  integrated 
internal  energy, 


J  dE  -  J  Cv  dT  +  [T  (3P/3T)v  -  P]  dv 


The  Hugoniot  pressure  can  be  expressed  as, 


P  =  P  +  Lg  (x)/C  ]  [P  v  (1  -  1/x)  -  F  (x)]  [l  -  v  (1  -  l/x)  g  (x)/2C  ]"  (25) 

ri  o  V  o  0  o  v 


where  PH  is  the  Kugoniot  pressure  and  P^  is  the  isothermal  pressure  at  the  same  volume 
ratios  v  /v.  The  other  quantities  are  defined  below 


X  =  V  /V 

o 


g(x)  =  B  a  x 
o 


F(x)  -  B  v  [  (1  -  s  or  T  )  (x  -  1  )/s  (s  -  1)  -  (1  -  l/x)/s. 
co  o 


Equation  25  is  based  on  the  Mumaghan  Equation.  The  Birch  and  Tait  Equations  can  also  be 
employed  in  such  derivations,  but  the  Murnaghar  Equation  gives  the  simplest  expression  for 


The  shock  compression  data  for  several  vinyl  polymers  are  available''  ;  however,  except 
for  density,  other  characteristics  of  the  polymers  used  in  those  experiments  and  the  initial 
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conditions  in  the  experiments  are  not  known.  In  order  to  use  Equation  -5,  one  must  select 
constants  in  the  Murnaghan  equation  for  the  given  material  and  given  initial  temperature  of 
the  shock  experiment. 

/o/?  \ 

It  is  stated  ,  that  the  material  used  in  the  shock  compression  of  polychlorotrifluoroethylene 

(Kel-F)  has  an  initial  density  p  =  2. 1,  but  the  initial  temperature  of  the  experiment  is  not 

°  3  (5) 

given.  For  this  density  the  specific  volume  is  0.  4762  cm  /g.  The  data  of  Weir  and 

(23) 

Hoffman  and  Weeks'  '  ’  indicate  that  this  material  has  a  degree  of  crystallinity  of  about  26 
percent.  Using  the  data  in  Reference  23  and  Reference  5,  degree  of  crystallinity  is 
calculated  for  a  series  of  specific  volumes.  With  this  information,  the  parameters  in  the 
Murnaghan  Equation  as  functions  of  crystallinity  can  be  plotted  and  extrapolated  to  the 
degree  of  crystallinity  of  the  test  sample.  Thus,  the  Murnaghan  parameters  s  and  are 
oo:cined  as  7.4  and  10,761  atmospheres,  respectively. 

In  the  case  of  plexiglas  (it  is  assumed  to  be  polymethylmethacrylate),  the  p^  given  was  1.  IS 

g/c.  c. ,  corresponding  to  a  specific  volume  of  0.  8475  cm3/  g.  By  use  oi  the  specific  volume 

\ 

as  a  function  of  temperature  at  one  atmosphere  reported  by  Hellwege  et  al.  ,  it  is  estimated 

that  this  p  corresponds  to  a  temperature  of  42.  5°C  (108.  5°F).  From  Figure  1  the 

°  2 
parameters  s  and  B^  for  the  Murnaghan  Equation  are  found  to  be  8.  70  and  35,550  Kg/cirT, 

respectively. 

In  Reference  36,  it  is  stated  that  the  polyethylene  used  in  a  shock  compression  test  has  a 
density  of  0.  92  g/c.  c.  T'hree  materials  of  the  same  density  (0.  92  g/c.  c. )  and  differing 
molecular  weight,  between  14,  000  and  33,  000  (determined  by  viscosity),  was  used  by 
Bridgman^10^  for  static  compression  to  40,  000  Kg/on,2.  The  branched  polyethylene  used  by 
Weir^3  has  a  density  of  0.  918  g  c.  c.  at  20°C.  A  high  pressure  polyethylene  with  a  density 
of  0.  918  g/c.  c.  at  20°C  was  also  used  by  Hellwege  et  al.  However,  it  is  evident  from 
Figure  16  that  the  pressure  versus  v  /v  curves  for  all  static  compression  lie  above  that 
for  the  shock  compression  test.  For  a  given  material,  however,  this  situation  is  thermo¬ 
dynamically  impossible.  This  seems  to  indicate  that  a  material  of  much  higher  compressibility 
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has  been  used  in  the  hock  experiment,  or  the  ambient  temperature  at  the  site  of  experiment 
is  much  higher  than  the  temperature  at  which  the  density  was  measured. 
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7.  DISCUSSION  AND  CONCLUSIONS 


The  results  of  our  recent  interest  in  the  equations  of  state  of  high  polymers  has  been  reported. 
Adoption  of  the  semi-empirical  equations  of  state  for  solids  for  use  with  polymers  seems  to 
have  a  wide  range  cf  applicability.  The  results  presented  in  this  paper  indicate  that  they 
are  capable  of  predicting  the  isothermal  volumetric  data  with  good  accuracy  to  very  high 
pressures.  They  give  smooth  and  reasonable  values  for  the  first  and  second  derivatives. 

The  compressibility  of  polymers  calculated  from  these  equations  are  quite  close  to  the 
values  obtained  from  graphically  determined  slopes  and  the  values  calculated  from  empirical 
polynomial  equations. 

The  ability  of  these  equations  to  predict  static  data  (Table  5)  at  high  pressures  much 

2 

beyond  the  original  data  range  is  quite  impressive.  The  data  calculated  to  40,000  Kg/cm 

in  Table  5  were  based  on  constants  which  were  obtained  from  the  data  of  Hellwege  et  al. 

2 

up  to  only  2000  Kg/cm  . 

By  use  of  the  theory  of  Slater^^  or  the  generalized  theory  by  Al'tshuler^^,  the 
Gruneisen  parameters  can  be  obtained  fi'om  these  equations.  The  Mumaghan  Equation 
yields  Gruneisen  parameters  dependent  only  on  temperature;  while  the  Tait  and  generalized 
Birch  equations  both  predict  Gruneisen  parameters  dependent  on  both  volume  and  temperature. 
The  values  thus  predicted  seem  to  possess  the  correct  theoretical  magnitude  at  normal 
temperatures. 

When  combined  with  the  Kennedy  fusion  lav/,  these  equations  predict  the  melting  temper¬ 
atures  of  polyethylene  (Table  9)  to  within  10  percent  at  pressures  mere  ihan  fifteen  times 
the  original  pressure  range. 

(5 


Under  shock  compression,  the  temperature  is  very  high  but  the  polymers  remain  solid 
because  of  the  high  melting  point  under  pressure.  Normalization  of  the  temperature  scale 
by  the  melting  temperature  may  extend  the  range  of  applicability  of  the  parameters  in  the 
equations  of  state  to  the  solid  in  this  high  temperature,  high  pressure  region. 

The  range  of  applicability  and  the  accuracy  of  these  equations  of  state  may  be  enhanced  if 

they  are  fitted  to  accurately  measured  compressibility  or  bulk  modulus  data  rather  than 

specific  volume  data.  As  proposed  by  Anderson  v  ’  '  r  easurements  of  ultrasonic  speed 

in  materials  under  high  pressure  is  most  suitable  for  obtaining  accurate  compressibility 

[  -  (bln  v/5  P)T1  or  bulk  modulus  [  -  (3  P/5  In  v)T]data.  The  compressibility  of  polymers 

has  values  in  the  order  of  10  3  which  cannot  be  obtained  accurately  from  direct  volumetric 

measurements.  Since  the  compressibility  is  the  first  derivative  of  (In  v)  with  respect  to  P, 

accurate  values  of  specific  volume  can  be  obtained  by  integrating  the  compressibility  data. 

From  such  data,  accurate  vakv  '  of  thermal  expansion  coefficient  of  polymers,  which  is 
-4 

in  ti.e  order  of  10  ,  can  alsc  a  obtained. 

/gg\ 

Guenther  et  al. v  '  have  measured  the  ultrasonic  speed  in  polymers  as  a  function  of  tem- 

(39) 

perature  under  atmospheric  pressure.  Zosaiv  " '  measured  the  ultrasonic  speed  in 
polyvinylchloride  at  1  atmosphere  and  1000  atmospheres  from  20°C  to  MO  C.  If  accurate 
values  of  are  available,  accurate  values  oi  compressibility  under  atmospheric  pressure 
can  be  obtained  as  a  function  of  temperature  irom  Reference  (38).  A  mere  detailed 
description  and  specification  of  the  polymers  and  the  imtial  conditions  of  the  experiment 
other  than  just  a  statement  of  the  density  \  ould  be  valuable  for  theoretical  work. 

The  success  oi  ’pplying  to  polymers  P\e  semi  °mn-:-  '.cal  equations  based  on  finite  strain 
theory,  may  be  attributable  to  the  increase  in  dej  re**  of  crystallinity  in  polymers  under 
high  pressure.  Although  the  Tait  equation  has  no  to  oretical  basis,  it  nevertheless  gives  a 
bulk  modulus  very  close  to  the  values  calculated  from  the  nther  two  equations.  The  Tait 
Equation  constant  A  was  found  to  vary  with  temperature.  This  variation  becomes  large 
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when  the  melting  point  is  approached  (Figure  3).  It  is  speculated  that  the  Murnaghan 
parameter,  s,  and  the  parameter  a  in  the  generalized  Birch  equation  may  also  be  nearly 

it 

uni\ersal  constants  for  a  wide  range  of  materials. 

In  shock  compression,  temperature  rises  with  pressure.  Degree  of  crystallinity  of  polymers 
increases  with  pressure^  ’  but  decreases  with  increasing  temperature'  .  At  low 
temperatures,  the  degree  of  crystallinity  decreases  very  slowly  with  increasing  temper¬ 
ature.  As  the  melting  *cmpe~ature  is  approached,  the  degree  of  crystallinity  drops 
rapidly  to  zero  at  the  melting  point.  In  this  region,  the  rate  of  increasing  crystallinity 
produced  by  pressure  may  become  insignificant  in  comparison  to  the  rapid  decrease  of 
crystallinity  produced  by  temperature.  This  situation  may  account  for  the  decrease  in 
slope  of  the  Hugoniots  at  the  highest  shock  pressures. 
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NOMENCLATURE 


A 


A' 


a 


B 

o 


B 


B' 


c 


C 


P 


C 


v 


E 


E 


o 


F(x) 


Tait  Equation  Parameter,  Dimensionless. 

Defined  in  Equation  (4). 

Constant  in  the  Simor  Equation,  Equation  (18). 

Constant  in  the  generalized  Eirch  Equation,  Equation  (3). 

Constant  in  the  Mumaghan  Equation,  Equation  (2),  atm., 
or  Kg/cm^. 

Bulk  modulus  defined  as  -v  Op/ov)  T;  Constant  in  the 
Tait  Equation. 

Defined  in  Equation  <4). 

Constant  in  the  Simon  Equation,  Equation  (18). 

Specific  heat  at  constant  pressure,  cal/g,  °K. 

Specific  heat  at  constant  volume,  cal/g,  °K. 

Internal  energy. 

Internal  energy  at  a  reference  state. 

Defined  in  Equation  (25). 
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Defined  in  Equation  (25). 


g(x) 

Ko 


k 


P 


H 


P 

P 


o 


s 

T 

t 


Tm,  tm,  T°,  t° 
m  m 


T 


o 


v 


V 

0 


X 


Constant  in  the  generalized  Birch  Equation,  Equation  (3), 
atm. ,  or  Kg/cm^. 

Constant  in  the  Kennedy  fusion  equation,  Equation  (20). 

Hugoniot  pressure,  atm..  Equation  (25). 

Pressure  of  a  reference  state,  a'm.,  Equation  (9). 

Pressure,  atm. 

Melting  pressure,  and  reference  melting  pressure,  usually  1  atm. 

Constant  in  the  Mumaghan  Equation,  dimensionless. 

o 

Temperature,  K. 

Temperature,  °C. 

Melting  temperature,  and  reference  melting  temperature. 

o 

Reference  temperature,  K. 

3 

Specific  volume,  cm'  /g,  under  pressure  P. 

Initial  specific  volume,  at  1  atm. 

Defined  as  (vQ/v) . 


z 


Defined  ir.  Equation  (12). 
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GREEK  SYMBOLS 


a 

P 

T 

v 

P 


Thermal  expansion  coefficient,  defined  as 
(Mn  vAT)  (0K)_1 

Compressibility,  defined  as  -v  ^(^v/^p)^,  atm  \ 

Griineisen  parameter  defined  in  Equation  (8)  and 
Equation  (9). 

Frequency  of  harmonic  oscillators,  cps. 


Density  at  a  reference  temperature,  T0, 
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